Charge and spin inhomogeneous phases in the Ferromagnetic Kondo Lattice Model 
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We study numerically the one-dimensional ferromagnetic Kondo lattice. This model is widely 
used to describe nickel and manganese perovskites. Due to the competition between double and 
super-exchange, we find a region where the formation of magnetic polarons induces a charge-ordered 
state. This ordering is present even in the absence of any inter-site Coulomb repulsion. There is 
an insulating gap associated to the charge structure formation. We also study the insulator-metal 
transition induced by a magnetic field which removes simultaneously both charge and spin ordering. 



I. INTRODUCTION 

A new paradigm in condensed matter springs up as 
a consequence of the growing body of experimental evi- 
dence suggesting different types of inhomogenous states 
in strongly correlated systems with competitive interac- 
tions. The simultaneous appearence of charge and spin 
superstructures seems to be present in completely differ- 
ent systems. For example, charge, lattice and spin mod- 
ulations with a doping-dependent wave vector have been 
observed in La2- a; Sr a ;Ni04 by means of neutron scatter- 
ing [1] and electron diffraction [2]. Stripe formation, to- 
gether with incommensurate spin fluctuations in high Tc 
superconductors can also be regarded as a manifestation 
of similar phenomena [3] . The charge and spin ordering 
found in many of the doped manganese perovskites fall 
in the same category. Experiments have revealed charge 
ordering (CO) at half filling in Ndo.5Sr .5Mn0 3 [4] and 
similar compounds, such as Pro.sCao.sMnOa, [5-8]. More 
recent interest has focused on electron doped charge or- 
dered manganites. The first observation of charge, or- 
bital, and magnetic ordering in electron doped mangan- 
ites dates back to the work of Wollan and Koehler [9] on 
Lai-^CazMn O3. Several works confirmed the observa- 
tion of this behavior in other compounds [10-12]. CO has 
also been found for other dopings as in Bii-^Ca^MnOa 
[10] and in Lai-^Ca^MnOa (doped with P r ) [12] for 
x>\ (few electron region). 

Manganese oxides have been intensively studied dur- 
ing the last years due to their collosal magnetoresistance 
(CMR) effect, i.e. a large change in the resistivity ac- 
compained by a transition from an insulating antiferro- 
magnetic (AF) or paramagnetic state to a ferromagnetic 
(FM) metal. In the above mentioned cases, the AF state 
also exhibits charge ordering, and the insulating gap is 



associated to the formation of the charge structure. 

Numerical studies for several dopings [13] have in- 
cluded the effect of nearest neighbor density-density re- 
pulsion V and obtained a very rich phase diagram. The 
diagram includes phase separation for either extreme 
dopings (1 ~ or 1 ~ 1), and CO for concentration 
x ~ 0.5. Using Monte Carlo techniques, Yunoki et. al. 
[14] studied the half filled case (x = 0.5) including both 
e g orbitals. They found that the CO phase is stabilized 
by Jahn- Teller phonons. However, as it is mentioned in 
their work [14], the z-axis stacking of charge and the ex- 
istence of bistripes at x > 0.5 [11] are both penalized by 
a large Coulomb interaction. This observation indicates 
that the Coulomb repulsion V is smaller than the value 
which is required to drive a charge ordered state. 

Unlike non-magnetic materials, where the Peierls in- 
stability or large inter-site Coulomb interactions are re- 
quired to get CO, we show here that charge density waves 
can result in manganites from the formation of mag- 
netic superstructures arising from the presence of com- 
peting interactions. A very rich family of inhomogeneous 
spin and charge structures [15] can be obtained with this 
mechanism, by changing the carrier concentration or the 
ratio between the competing interactions. 

In order to illustrate our new mechanism, we will study 
here a simplified model where the competing forces are 
represented by the double and super-exchange interac- 
tions. The Ferromagnetic Kondo Lattice Model (FKLM), 
proposed to describe the manganites, was studied orig- 
inally by de Gennes [16]. One of the most intriguing 
questions related to this model is: what arc the emerging 
phases when the competing interactions are energetically 
similar? De Gennes proposed canting of two interpene- 
trating lattices as the compromise solution for this com- 
petition. This concept was used subsequently in recent 



1 



analytical approaches [17,18]. In addition, phase separa- 
tion has also been considered as a possible solution to this 
competition [19,20]. In contrast, in the present paper we 
find spin phases which cannot be described in terms of 
two interpenetrating sublattices nor do they correspond 
to phase separation. 

As a result of the competition between the double 
exchange (DE) mechanism which delocalizes the hole, 
and the super-exchange (SE) between local spins, dif- 
ferent phases may appear such as ferromagnetism in one 
extreme, antiferromagnetism in the other, and doping- 
dependent modulated charge and spin order in between. 
For x = 1/2, we do not find charge ordering for V = 0. 
However for concentrations which are away from half fill- 
ing, we obtain charge modulation even in the absence of 
Coulomb interactions. This new state can be described 
as an ordered phase of FM islands where carriers are lo- 
calized. 

Since the charge modulation is driven by the AF or- 
dering between the islands, the transition to a FM state 
destroys the charge ordering, closing the insulating gap. 
This transition can be induced by applying an external 
magnetic field. If the energy difference between the AF- 
charge ordered phase and the FM state is small, the 
transition is induced by a magnetic field which can be 
generated in the laboratory. In this respect, this new 
mechanism gives a natural way to understand how the 
interplay between charge and spin degrees of freedom can 
induce a collosal magnetoresistive effect. The complete 
understanding of this phenomena in manganites requires 
a more detailed study, exploring regions of concentration 
where there is no charge ordering since the concentration 
of particles is not commensurated with the lattice. 

To study the insulator-metal transition we have com- 
puted the Drude weight by exact diagonalization of finite 
size systems. The phase diagram has been computed us- 
ing the Density Matrix Rcnormalization Group (DMRG) 
[21,22]. The energy of the different phases has been com- 
puted exactly for the case where the localized spins are 
classical. In this way we compare the phase diagrams 
corrcponding to the classical and quantum limits of the 
spin variables. 

In the next section we introduce the model. In Sec- 
tion III, we study the phase diagram making a detailed 
description of its different phases and the way they have 
been obtained. In Section IV we study the charge exci- 
tations by calculating the Drude weight, the charge gap, 
and the effect of a magnetic field on the metal-insulator 
transition. Section V is devoted to the conclusions. 

II. THE MODEL 

In the manganese oxide compounds, the Mn 4+ ions 
consist of three localized t% g orbitals that can be repre- 
sented by an S — 3/2 spin due to the large Hund inter- 



action Jh- The Mn 3+ ions have an additional electron in 
the doubly degenerated e g orbital (where the degeneracy 
is lifted due to Jahn- Teller interactions). To describe this 
system we use the so called Ferromagnetic Kondo Lat- 
tice Model (FKLM) [23-25]. This Hamiltonian includes 
a lattice of localized spins which are coupled via Jh to 
itinerant e g electrons: 

<ij>cr i <ij> 

-J H J2Sm + Kj2SiSj (!) 

i <ij> 

where rii a — <^ ia Ci a and rii = njf+n^. The first term rep- 
resents the e g -electron hopping between nearest-neighbor 
Mn ions at sites i and j. The second and third terms are 
the on-site and nearest-neighbor Coulomb repulsions for 
the same orbitals. Jh is the Hund's coupling between 
localized Si, and itinerant <7j spins. We will denote the 
concentration of holes by x, i.e. (1 — x) is the concentra- 
tion of itinerant electrons. Hund's interaction together 
with the hopping term give rise to the DE interaction 
which favors FM ordering of the local spins. K is the an- 
tiferromagnetic SE interaction between local spins (t^g 
electrons). This term stabilizes the AF phase for x = 1, 
and competes with the DE term for intermediate fillings. 
For completeness we add the effect of a nearest-neighbor 
repulsion V when considering x = 0.5. The localized 
spins are taken as S = 1/2, instead of S = 3/2, to sim- 
plify the calculations. The model is studied in one dimen- 
sion (ID), since the basic features of the phase diagram 
for the 3D case appear in one and two dimensions as well 
[24,26]. 

Experiments indicate that Jh S> t. For this reason 
we fix the values of Jh = 20i and U = lOt. Then, the 
only free parameters are V and K. For the large value 
of Jh used in this paper, the results are not sensitive to 
U since double occupancy is suppressed by Jh- For this 
reason, we propose to use a t-J model for the conduction 
electrons hence reducing the Hilbert space by eliminating 
the doubly occupied states: 

<ij><y <ij> 

+V J2 mnj - J H J2 SitT-i + K SiSj 

<ij> i <ij> 

where c\ a = c\ a {\ - n lS ) and c ia = (1 - n m )c\ a , are 
the creation and anihilation operators of constrained 
fcrmions (no more than one electron per site). Strictly 
speaking, the exchange between two e g electrons depends 
on the configuration of the localized spins, involving 
a more complicated form of the reduced Hamiltonian. 
However, we are concerned with the low-energy prop- 
erties so that simply assuming the above form (2) with 
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J = u+ Jh is sufficient. From the expression for J, we see 
that the condition for the validity of H as a low energy 
effective theory for H is t <C U + Jh- We have checked 
on small systems that H gives similar results as the full 
Hamiltonian H. We will use H to calculate the Drudc 
peak and the charge gap (see Sec. IV). 

These models are solved using exact diagonalization 
(Lanczos) and DMRG techniques in open systems. When 
we used DMRG, we took special care in the growing pro- 
cedure to avoid lattice frustration. For instance, we con- 
sidered different right and left blocks. We obtained max- 
imum discarded weights of the order of 10~ 5 using the 
finite-size algorithm and keeping 150 states. To charac- 
terize the different phases we calculated static correlation 
functions for system sizes up to L = 40 sites. This is de- 
scribed in the following section. 

III. PHASE DIAGRAM 

The charge and spin correlation functions are com- 
puted for the complete Hamiltonian H (Eq. (1)), using 
DMRG and open boundary conditions. These correlation 
functions are the spin and charge structure factors: 

S{q) = \ j Y,e i(Ri - Ri)q (S&i) (3) 

N(q) = e^- R ^ ((m n)(nj n)) , (4) 

where n = 1 — x is the average number of electrons per 
site, L is the chain length, and Si is the localized spin 
operator at site i. Both correlation functions are used to 
build the phase diagram of H. We describe below the 
different phases obtained in this way. 

A. Ferromagnetic phase (FM) 

When K <t the local spins are aligned in a FM state 
(Figs. 1, 3a). At the cost of the AF interaction, the 
system gains kinetic energy when neighboring spins are 
parallel. In this situation, the itinerant electrons can be 
mapped into spinless fermions because the large Hund's 
interaction and the FM ordering of the localized spins 
force the band electrons to be fully polarized. 

The charge correlation function for the spinless model 
can be calculated analytically. It has the shape of a trun- 
cated pyramid. The truncation begins at q — 2kp, where 
for the spinless case kp — tt(1 — x). This is shown in 
Fig. 2 for x — 1/2. The oscillations are due to the 
open boundaries conditions. For this charge or- 

dered phase sets in when V reaches the critical value 



V c = 2t. There is a region above this critical value where 
the spins remain ferromagnetically ordered. However, 
larger values of V inhibit the DE mechanism by localiz- 
ing the charges, and other spin orderings arise (Fig. 1). 
For V = 0, our numerical calculations for the fully quan- 
tum case (S=l/2) lead to an upper boundary of the FM 
region at K=0.2t. 

B. Island phase (ILP) 

When the localized spins are considered as classical de- 
grees of freedom, canted or spiral states have been pro- 
posed [16] for the transition between the FM (K/t = 0) 
and the AF states (fx). These phases are obtained by 
variational methods which exclude inhomogeneous states 
as candidates to be the ground state. A more complete 
calculation with classical spins [27] shows that the 2fcp 
instability of itinerant electrons in one dimension leads 
to the type of spin and charge ordering that we find with 
our numerical calculations with quantum spins. Indeed 
our DMRG results confirm and extend this picture. For 
the intermediate K/t regime, the spin and charge struc- 
ture factors peak at q = fcp and q = 2kp respectively (Uf 
corresponds to the spinless case). 

The following are the results for some particular fill- 
ings: 

a) x = 2/3. In this case the phase roughly consists of 
FM islands of 3 spins which are aligned antiferromagnet- 
ically (TTTIIITTTIII)- The spin structure factor for the 
fully quantum case is shown in Fig. 3b. Two well defined 
peaks are observed at q = ±7r/3. The nearly logarithmic 
increase of these peaks with the system size L (see Fig. 
4b) indicates a power-law decay with exponent one for 
the spin-spin correlation as a function of the relative dis- 
tance. However the system sizes considered here do not 
exclude other possible values for this exponent. In any 
case, the linear increase of the peak in a logarithmic scale 
indicates that the correlations are quasi long-ranged. The 
absence of real long-range order is a consequence of the 
one-dimensional character of the model and the SU (2) 
continuous symmetry. The same behavior is observed 
for the other spin structures described below. 

In the classical limit there is no effective hopping be- 
tween spins which are antiferromagnetically aligned since 
tdass — tcos(0/2) (where 9 is the relative angle between 
the spins). Therefore it is clear that the charge is then 
localized in each island when the spins are classical. For 
the quantum case the charges are still localized because 
the effective hopping between spins in the same island is 
still larger than the hopping between spins in different 
islands. The last observation plus the fact that there is 
one electron per island allows to understand the opening 
of a gap at q = hp. 

The localization of the charge is confirmed by the non- 
uniform distribution shown in Fig. 5c. The translation 
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symmetry of H is spontaneously broken in the thermo- 
dynamic limit for the charge ordered state. Even though 
the symmetry is only broken for the infinite-size chain, 
the DMRG process selects one of the broken symmetry 
states also for finite chains. The charge occupations per 
site (0.45 in the borders and 0.27 in the center) are very 
close to the values expected for the classical limit (0.25 
in the borders and 0.5 in the center for a three sites open 
chain) . Consistently with this picture, the charge-charge 
structure factor (Fig. 2b) shows remarkable peaks for 
q = ±27r/3. Fig. 4a shows that the intensity of these 
peaks scales linearly with L. This is evidence of long- 
range charge ordering. This behavior is also found for 
other commensurated fillings like x = 1/3. 

b) x = 1/3. The classical image of the intermediate 
K/t phase for this filling is: TT4TTITTITT4- Here we find 
two-site FM islands separated by an antialigned spin (see 
Fig. 5a). S(q) shows peaks at q = ±27r/3 (see Fig. 3b). 
From the real space spin-spin correlation function we can 
rule out the uniform spiral phase since the 'islands' are 
clearly distinguishable (this is also valid for x = 2/3). In 
a classical image, the charge distribution should present 
one electron in the down spin and another distributed be- 
tween the two parallel spins. Our numerical results show 
a density of approximately 0.6 electrons on the parallel 
up spins (FM spins) and 0.8 on the anti-aligned spins. In 
this case, the inclusion of a hopping between the up and 
down spins is essential to describe the numerical results. 
N(q) shows the expected peaks at q = ±2tt/3 (Fig. 2b). 
The peak near zero is due to finite-size effects. 

One interesting point is that for classical spins there 
is electron-hole symmetry; i.e. x and (1 — x) should ex- 
hibit the same phase for a given set of parameters. The 
absence of that symmetry for quantum spins is reflected 
in our results (the ILP for x = 1/3 is different from the 
one corresponding to x = 2/3). 

c) x = 1/2. In this case the FM islands contain two 
spins and, like in the previous cases, the islands are an- 
tifcrromagnetically ordered: (Ttilttii)- S(q) shows a 
clear peak at q = ±7r/2 as it is expected for a unit 
magnetic cell containing four spins (see Fig. 3). The 
same structure has been obtained with the Monte Carlo 
method for classical spins [25]. However, this state has 
been interpreted as a spiral phase. Other classical Monte 
Carlo calculations [28] confirm the peak at q — tt/2 but 
instead of a spiral phase, they show the island state that 
we obtain with quantum spins. The effective hopping al- 
ternates between two values from one bond to the next 
one in this particular island phase being larger between 
the parallel spins. This alternation of the hopping does 
not induce any charge inhomogeneity, i.e. if V = there 
is no charge density wave for this concentration (see Figs. 
4, 5 and 6). However N(q) differs from the FM case. 
This change is related to the presence of bond order- 
ing (dimerization) . The charges are 'localized' in bonds 
due to the particular spin structure generated by the two 



competitive interactions. It is well known that this type 
of ordering also opens a charge gap which gives rise to 
an insulating behavior. We will come back to this in Sec. 

IV. The shape of N(q) is closer to (1 — cosg)/4 (charge 
structure factor of a completely dimerized state) as a con- 
sequence of the enhancement of the charge correlations 
within the islands (bonds). 

The above results are for V = 0. When increasing the 
Coulomb interaction V, a CO phase appears for V ~ 2t 
like in the spinless model. For V > 2t, the usual CO 
phase with a charge structure factor peaked at q = ir 
is obtained. S(q) still has a peak at q — tt/2 (Figs. 
1, 7). Larger values of V inhibit the double exchange 
mechanism since the charge fluctuation are strongly su- 
pressed. The magnetic correlations change then from 
q = tt/2 to a ferrimagnetic (FIM) CO phase. This phase 
can be schematically represented by 1U1U1U, where the 
total spin per site is alternating between S = l(fr) and 
S = 1/2 (J.). The S = 1 spin at site i corresponds to 
the local triplet state forced by Jh between the itinerant 
electron and the localized spin at the same site. 

The same calculations for x = 4/5 and 3/4 show spin 
and charge correlation functions consistent with the for- 
mation of five and four FM-spin islands respectively. 
These structures could be regarded as a crystallization 
of the magnetic polarons described in references [29] for 
the dilute limit. It is interesting to note that charge or- 
dering is induced by spin ordering and vice versa. The 
formation of these superstructures is then a consequence 
of the interplay between charge and spin degrees of free- 
dom. 

C. Antiferromagnetic phase (AF) 

As K is increased, the AF ordering prevails over the 
DE mechanism for any concentration. For instance, at 
half-filling the transition to the AF phase occurs for 
K/t ~ 0.7 (for x ^ 1/2 the transition occurs for smaller 
values of K/t). S(q) shows a strong peak for q = ir and 
the charge is uniformly distributed along the sample, i.e. 
N(q) is broader. As mentioned above, a FIM charge 
ordered state is obtained for half filling by increasing 

V. (Figs. 1,7). The charges remain in one sub-lattice 
forming a chain with alternating S=l and S=l/2 spins 
which are antiferromagnetically aligned. 

The validity of the picture described above can be 
checked in the classical limit for the localized spins. We 
calculated the contributions of DE and SE to the en- 
ergy of the different possible phases for each concentra- 
tion x. To this end, we assumed that the hopping at 
the AF bonds is tAF — t/\/2 (the conclusions are the 
same if we take tAF = 0). This value of tAF gives a 
better agreement for the charge distribution in the ILP 
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states and with the critical values of K/t separating dif- 
ferent phases. For instance, for i = 1/2 we obtained 
E FM = -2t/Tr + K/4-J H /8 , = -.564761*- J H /8, 
and Eaf — ~V2t/ir — K/A — Jh/8. In this way, the se- 
quence of stable spin phases as a function of increasing 
K/t, FM^ 7r/2 — »AF, coincides with our numerical re- 
sults for quantum spins (see Figs. 1 and 8b). This simple 
calculation clearly shows how the FM state weakens by 
increasing K/t and gives rise to more complicated struc- 
tures when both interactions are comparable. In this 
classical picture, a canted AF phase has a slightly lower 
energy than the pure AF phase. However the 7r/2 is- 
land phase remains as the most stable state over a finite 
interval in K/t. The same procedure leads to similar 
conclusions for x = 2/3 and 1/3 (see Figs. 8a and c). 

D. Magnetic mechanism for pairing 

As K is further increased (K > t and V = 0), we find 
a peak at q — n/A in the charge structure factor(Fig. 6). 
To understand the origin of this peak we will consider the 
effective Hamiltonian which emerges by taking the limit 
of infinite Jh- In this limit, the original Hamiltonian H 
(Eq. (1)) can be projected into the subspace having fi- 
nite energy, i.e., each carrier represented by an effective 
spin S = 1 and each hole represented by a spin 5 = 1/2. 
Double occupation is projected out for infinite Jh . The 
projected Hamiltonian can be written as: 

He ff = t E Pi,j(2 n irij -m- rij)(SiSj + ^) 

<i,j> 

+ KJ2 [0--ntKl-n j )(S i S j -\) 

<i,j> 

+ ^niUj -rii- rij)(SiSj - ^) 

+ ^.(SiS,- - 1)] (5) 

where i,j are nearest neighbors and rii is the occupation 
number at site i. rij = if there is a spin S — 1/2 at 
site i (hole) and rii = 1 if there is a spin 5 = 1 (one 
conduction electron). Pij is a permutation operator be- 
tween the spins at sites i and j. Note that this permu- 
tation is operating only between spins 1 and 1/2 due to 
the factor (2rnnj — rii ~ n j)- This is because the hop- 
ping process is only possible if the occupation numbers 
at sites i and j are different. The last three terms of Eq. 
(5) are derived from the original AF interaction K be- 
tween localized spins. This interaction remains the same 
between two spins 5 = 1/2, but the situation is different 
between a spin 5=1 and a spin 5 = 1/2, or between 
two spins 5=1. In the last two cases the effective AF 
interaction is reduced due to the fact that the localized 
spin 5=1/2 which is part of a spin 5=1 cannot form 



a singlet state with its nearest neighbor spin. While the 
minimum AF energy for the singlet state between two 
spins 5 = 1/2 (2 holes) is —K, the corresponding en- 
ergy for the doublet between a spin 5 = 1 and a spin 
5=1/2 (one particle and one hole) is — The singlet 
between two spins 5=1 (two particles) has an AF energy 
— The renormalization of the AF interaction gives 
rise to an attractive interaction for the charge degrees of 
freedom. The magnitude of this interaction is given by 

^e// = E2part+E2holes~-2Ei part — 3 K+— = — — . 

It is well known that a nearest-neighbor attraction gives 
rise to three different phases depending on the ratio be- 
tween the attractive interaction V e ff and the effective 
hopping (~ t/2). If this ratio is smaller than some criti- 
cal value (~ \/2) the system is a Luttinger liquid, and the 
metallic correlations are the dominating ones [30] . Above 
that value the system remains as a Luttinger liquid but 
the superconducting correlations become dominant. It is 
precisely in this region where the formation of pairs gives 
rise to a peak at q = 7r/4 in the charge structure factor 
for half filling. If the attractive interaction V e // is larger 
than the double of the effective hopping (V e ff > t), the 
system segregates into two different phases. One phase is 
rich in spins 5 = 1/2 (holes) and the other one is rich in 
spins 5=1 (electrons). A first indication of this behav- 
ior is presented in [23] , where binding and phase segrega- 
tion is obtained in a similar model. However the binding 
energy is overestimated because the authors have only 
included the AF interaction between spins 5 = 1/2, i.e., 
they did not renormalized the AF effective interaction. 

IV. METALLIC AND INSULATING PHASES 

The diversity of phases found for the FKLM sug- 
gests different transport properties. To determine the 
metal or insulator character of the phases we have calcu- 
lated the charge gap and the Drude weight (V = 0) for 
x = 1/3, 1/2, and 2/3. 

The charge gap for systems having N particles is de- 
fined as: 

A c = Ejy + i + Ejy-i — 2En. (6) 

En is the ground state energy of the considered system. 
En + i and Em~i are the ground states for N + 1 and 
N — 1 particles respectively. To get the charge gap in the 
thermodynamic limit we did a 1/L extrapolation of the 
gap A C (L) for different sizes. 

The Drude weight D(L) is defined as: 

where <f> can be interpreted as the magnetic flux inside a 
closed chain. D(L) was calculated by exact diagonaliza- 
tion in closed chains having 8 and 12 sites. The ground 



5 



state used to compute D(L) is the one which minimizes 
the energy as a function of the flux <fr. To reduce the 
Hilbert space, we did this calculation with the Hamil- 
tonian (2). The results were tested by doing the same 
calculation with the full Hamiltonian (1) in a half filled 
chain with 8-sites. We did not find any substantial vari- 
ation. 

In the FM phase, the charge gap is zero as it is shown 
in Fig. 9. Consistently with the cancellation of the gap, 
there is a finite Drudc weight which coincides with the 
one obtained for a spinless model (see Fig. 10). When 
K/t is large enough to get into the ILP, D(L) decreases 
abruptly and approaches zero. This is indicative of a 
metal-insulator transition. These results are not con- 
clusive because we cannot solve large enough chains to 
make an accurate extrapolation to the thermodynamic 
limit. Another indication of the insulating character of 
the ILP state for x = 1/2 and 2/3 is the finite value of 
the charge gap. The gap for x = 1/2 is greater than the 
one for x = 2/3, as it is expected from a simple calcula- 
tion with classical spins. For x = 1/3, the gap calculated 
in the ILP is small (A c < 0.05i), but the precision does 
not allow us to obtain a reliable extrapolation to the 
thermodynamic limit. D increases for larger values of 
K, but it is unclear whether the AF state is insulating 
or metallic. The charge gap decreases and cancels when 
the AF phase is established (Fig. 9). The closure of the 
charge gap and the finite value of D suggest a metallic 
character for this phase. 

To investigate the possibility of an insulator-metal 
transition with magnetic field applied to the ILP phase, 
we have calculated the lowest energy state for different 
magnetizations (L < 40). For x = 1/2 we found that the 
total spin increases smoothly with field and the charge 
gap closes abruptly when the system reaches its maxi- 
mum magnetization (see Fig. 11). A similar behavior is 
seen in the ILP of x = 2/3. As it is expected, the field 
needed to induce the transition from the ILP to the FM 
phase decreases when K/t gets close to the critical value. 
This is shown in Fig. 8, where we compare the energies 
of the different phases. 

V. CONCLUSIONS 

In summary, we have presented numerical evidence for 
the existence of a new type of charge and spin ordering 
in the FKLM. This mechanism, induced by the compe- 
tition between double and super-exchange, is based on 
a striking interplay between charge and spin degrees of 
freedom. As we have shown in a previous work [29], in 
the dilute limit each carrier polarizes its surroundings in- 
ducing a FM distortion. The quasiparticles constituted 
by the electron plus the magnetic distortion are called 
FM polarons. The size of these polarons is governed by 



the ratio K/t. When the mean separation between carri- 
ers equals the size of the FM distortion associated to each 
polaron (FM islands), the SE interaction induces an AF 
ordering between polarons. This is of course a schematic 
picture since there are quantum magnetic fluctuations 
within and between the islands. 

While previous works reported the need of a Coulomb 
interaction V [13] or electron-phonon coupling A [31] to 
stabilize a charge-ordered phase, we show that it arises 
naturally just by considering the competition between 
double and super-exchange. 

A great variety of experiments have found si- 
multaneous charge and spin ordering in manganites. 
The extraordinary colossal magnetoresistance effect for 
Lao.sCao.sMnOs involves the abrupt destabilization of a 
CO-AF state by application of a magnetic field [32]. In- 
sulating charge-ordered and metallic FM regions coexist 
in (Lao.5Ndo.5)2/3Cai /3 Mn03 [33] and Pro. 7 Cao.3Mn0 3 
[34,35]. Both phenomena indicate that the CO phase is 
very close in energy to the FM state. We have shown 
here that a magnetic field can produce an abrupt metal- 
insulator transition by polarizing the local spins and 
favouring the double exchange mechanism. For these 
reasons, we think that the phenomenon presented here 
is relevant for the underlying physical mechanism stabi- 
lizing spin and charge structures in manganites. 

Our model calculations fail to include several effects 
that may play an important role in real systems like 
Jahn- Teller distortions and orbital degeneracy. However, 
the spirit of the present paper is to show that a mini- 
mal model, including the most relevant electron-electron 
interactions, is sufficient to explain the simultaneous ap- 
pearing of charge and spin ordering. This is a possibility 
which had been excluded by the spin canted phases pro- 
posed in previous works [16]. 
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FIG. 1. Phase diagram for x = 1/2. The circles indicate 
the points where calculations have been performed. CO, I 
and M stand for charge order, insulating and metallic phases 
respectively. FM, FIM, AF and n/2 stand for ferro-, ferri-, 
antiferromagnetic and q = -k/2 spin phases respectively. 




FIG. 2. Charge correlations for the different phases men- 
tioned in Fig. 1 (V = 0). 




FIG. 4. a) Charge and b) spin peak-height scalings in 
the island-type phases. The spin correlation is plotted on 
a semi- logarithmic scale. 
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FIG. 5. Charge mean values (bars, right axis) and spin 
real space correlations (squares, left axis) for the island type 
phases. The horizontal line represents the average charge. 




FIG. 3. Spin correlations for the different phases found 
(V — 0): a) Ferromagnetic phase (FM); b) Island-like Phase FIG. 6. On-site charge mean values for x=l/2 in the FM 

(ILP); c) Antiferromagnetic Phase (AF) (black), ILP (striped) and AF (empty) phases. 
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FIG. 7. Spin and charge correlations for half filling in the 
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FIG. 8. Classical energies for different fillings in the FKLM. 
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FIG. 9. Extrapolated charge gap in the thermodynamic 
limit. The lines are only a guide to the eye. 




FIG. 10. Drude weight for different fillings, J/t = 0.4, 
V = and L = 12 sites as a function of the super-exchange. 




FIG. If. Mean on-site magnetization S z and extrapolated 
charge gap (A c ) as a function of the magnetic field (H) in the 
ILP phase for x=l/2. 
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